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The coronavirus nucleocapsid (N) protein plays a multifunctional role in the virus life cycle, from regulation of replication and 
transcription and genome packaging to modulation of host cell processes. These functions are likely to be facilitated by interac- 
tions with host cell proteins. The potential interactome of the infectious bronchitis virus (IBV) N protein was mapped using sta- 
ble isotope labeling with amino acids in cell culture (SILAC) coupled to a green fluorescent protein-nanotrap pulldown method- 
ology and liquid chromatography-tandem mass spectrometry. The addition of the SILAC label allowed discrimination of 
proteins that were likely to specifically bind to the N protein over background binding. Overall, 142 cellular proteins were se- 
lected as potentially binding to the N protein, many as part of larger possible complexes. These included ribosomal proteins, 
nucleolar proteins, translation initiation factors, helicases, and hnRNPs. The association of selected cellular proteins with IBV N 
protein was confirmed by immunoblotting, cosedimentation, and confocal microscopy. Further, the localization of selected pro- 
teins in IBV-infected cells as well as their activity during virus infection was assessed by small interfering RNA-mediated deple- 


tion, demonstrating the functional importance of cellular proteins in the biology of IBV. This interactome not only confirms 
previous observations made with other coronavirus and IBV N proteins with both overexpressed proteins and infectious virus 
but also provides novel data that can be exploited to understand the interaction between the virus and the host cell. 


F or positive-strand RNA viruses, the viral genomic RNA-bind- 
ing protein (often referred to as a capsid or nucleocapsid [N] 
protein) plays a key role in virus biology, with functions in encap- 
sidation and assembly, and for many viruses contributes toward 
the regulation of viral RNA synthesis (e.g., see reference 1). These 
proteins may also interact with host cell proteins both to facilitate 
their activities related to the virus life cycle and possibly to disrupt 
host cell signaling pathways (e.g., see reference 2). Coronaviruses 
and arteriviruses are positive-strand RNA viruses with common 
genome replication and expression strategies that encode an N 
protein that has multiple functions in the virus life cycle (3-8). 
Many of these functions may be mediated by subcellular localiza- 
tion (9-11) and multiple interactions with host cell proteins. 
The coronavirus N proteins show high sequence homology 
within each genus but divergent homology between genera (12). 
This is illustrated by the conservation of functional regions within 
N proteins, such as RNA-binding motifs and membrane (M) pro- 
tein binding domains (13). However, the location of such sites in 
the primary structure of the protein differs between genera. The N 
protein is best known as a viral structural protein, whose function 
is to bind the RNA genome to form a helical ribonucleoprotein 
(RNP) in mature virions (14). The N-terminal region of N protein 
has been implicated in RNA binding for infection bronchitis virus 
(IBV) (15-17), severe acute respiratory syndrome coronavirus 
(SARS-CoV) (18), mouse hepatitis virus (MHV) (19), and human 
coronavirus OC43 (20), as have additional sites in region 2 of 
MHV (21) and region 3 of IBV (22). Structural data suggested 
conserved tertiary structures between these viruses which may be 
involved in RNA binding (23). Phosphorylation has been shown 


9486 _ jvi.asm.org 


Journal of Virology p. 9486-9500 


to play a role in the discrimination of viral and nonviral RNA in 
IBV N protein (15, 24), with glycogen synthase kinase 3 (GSK3) 
being the acting kinase for SARS-CoV and MHV N protein. IBV N 
protein is postulated to bind RNA via a lure-and-lock mechanism 
(15, 16). 

Aside from the structural functions, N protein has been impli- 
cated in several other processes. The SARS-CoV N protein has 
been shown to alter the activities of the AP-1, Akt, Erk, Jnk, and 
p38 mitogen-activated protein kinase signaling pathways as well 
inducing a block in the Go/G, phase of the cell cycle (25-27). The 
association of N protein with elements of the viral genome as well 
as cellular RNA-binding factors strongly suggests a role in the 
transcription of subgenomic mRNAs. N protein has been shown 
to be essential for the efficient rescue of a number of infectious 
clones of coronavirus genomes, including IBV (28). While repli- 
cation and transcription can occur in the absence of N protein, the 
efficiency of both processes requires N protein to be provided 
either in cis or in trans (4, 6, 29). 

N protein colocalizes with cytoplasmic replication complexes at 
early time points postinfection in both transmissible gastroenteritis 
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coronavirus and MHV and does so in a dynamic manner (30, 31), 
which agrees with the suggested involvement of the coronavirus N 
protein in virus replication (29). In a wide range of coronaviruses and 
arteriviruses, the N protein also localizes to the nucleolus (9, 32-34), 
a subnuclear structure whose key role is in ribosome biogenesis, but 
with additional roles in cell cycle regulation and the response to cell 
stress (35). Some variable degree of nuclear localization may also be 
observed, and many positive-strand RNA virus RNA-binding pro- 
teins have been shown to be located in the cytoplasm and nucleolus 
(for a review, see reference 36). 

The interaction of coronavirus N proteins with host cell factors 
have not been subjected to extensive study, and therefore, the 
interactome of a coronavirus N protein was investigated through 
the combination of stable isotope labeling with amino acids in cell 
culture (SILAC)-based mass spectrometry (MS) and green fluo- 
rescent protein (GFP)-nanotrap pulldown methodologies. IBV N 
protein was chosen because both its biophysical properties (15- 
17, 24, 37-39) and subcellular localization to the cytoplasm and 
nucleolus (40—44) have been extensively characterized, as have the 
changes in the host cell proteome as a result of infection with IBV 
(45, 46). 


MATERIALS AND METHODS 


Expression and harvesting of EGFP and EGFP-N in 293T cells. 293T 
cells were cultured in SILAC-Dulbecco’s modified Eagle medium 
(DMEM) lacking arginine (R) and lysine (K) supplemented with versions 
of these amino acids containing ROKO (for enhanced GFP [EGFP] expres- 
sion) and R6K4 (R = 'C-6, K = D-4) (for EGFP-IBV N expression). Cells 
were maintained in SILAC medium supplemented with 10% dialyzed fetal 
bovine serum (FBS) and 5% penicillin-streptomycin for at least 5 cell 
divisions to ensure >95% labeling of cells. For the passage of cells, 
trypsinization was avoided, with cell dissociation buffer (phosphate-buff- 
ered saline [PBS] based; Invitrogen) being used instead. Calcium phos- 
phate was used to transfect cells with plasmid DNA. Per 10-cm? dish, 500 
wl of 2X HEPES-buffered saline and 1.5 mM Na,HPO, (in which the pH 
was adjusted to 7) was added dropwise to a solution of 61 wl of 2 M CaCl, 
and 10 pg of DNA (diluted in nuclease-free water); this solution was then 
added to 293T cells. 

Lysis and immunoprecipitation of GFP-fusion proteins were per- 
formed using GFP-trap beads (Chromotek) that consist of a single-do- 
main anti-GFP antibody conjugated to an agarose bead matrix. Cell pel- 
lets were incubated for 30 min with 200 pl lysis buffer, 10 mM Tris-HCl, 
pH 7.5, 150 mM NaCl, 0.5 mM EDTA, 0.5% NP-40, and 1x EDTA-free 
protease inhibitor (Roche). The lysate was cleared by centrifugation and 
diluted 5-fold with dilution buffer comprising 10 mM Tris-HCl, pH 7.5, 
150 mM NaCl, 0.5 mM EDTA, and 1X EDTA-free protease inhibitor 
(Roche). The GFP-trap beads were equilibrated with ice-cold dilution 
buffer and then incubated with diluted cell lysate for 2 h at 4°C ona rotary 
mixer, followed by centrifugation at 2,700 * g for 2 min. The bead pellet 
was washed once with dilution buffer, followed by a single wash in buffer 
comprising 10 mM Tris-HCl, pH 7.5, 300 mM NaCl (concentrations of 
150 mM and 500 mM were also tested [see below]), 0.5 mM EDTA, and 
1X EDTA-free protease inhibitor (Roche). After centrifugation of the 
GFP-trap beads at 2,700 < g and removal of the wash buffer, the beads 
were resuspended in 2 SDS-sample buffer and boiled for 10 min to elute 
bound proteins. All stages of the processes were conducted on ice or 
at 4°C. 

LC-MS/MS. Protein samples generated by GFP-trap immunoprecipi- 
tations were separated by one-dimensional SDS-PAGE (4 to 12% bis-Tris 
Novex minigel; Invitrogen). The resulting separated proteins were cut 
from the gel in 10 slices and subjected to in-gel digestion with trypsin. 
Trypsin-digested peptides were separated using an Ultimate U3000 nano- 
flow liquid chromatography (LC) system (Dionex Corporation) consist- 
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ing of a solvent degasser, micro- and nanoflow pumps, a flow control 
module, a UV detector, and a thermostated autosampler. A sample vol- 
ume of 10 wl (comprising 2 wg) was loaded at a constant flow rate of 20 
l/min onto a PepMap Cj, trap column (0.3 mm by 5 mm; Dionex Cor- 
poration). After trap enrichment, peptides were eluted onto a PepMap 
C,g nanocolumn (75 ym by 15 cm; Dionex Corporation) with a linear 
gradient of 5 to 35% solvent B (90% acetonitrile with 0.1% formic acid) 
over 65 min at a constant flow rate of 300 nl/min. The high-pressure liquid 
chromatography system was coupled to an LTQ Orbitrap XL instrument 
(Thermo Fisher Scientific Inc.) via a nanoelectrospray ion source 
(Proxeon Biosystems). The spray voltage was set to 1.2 kV, and the tem- 
perature of the heated capillary was set to 200°C. Full-scan MS survey 
spectra (m/z 335 to 1,800) in profile mode were acquired in the Orbitrap 
instrument with a resolution of 60,000 after accumulation of 500,000 ions. 
The five most intense peptide ions from the preview scan in the Orbitrap 
instrument were fragmented by collision-induced dissociation (normal- 
ized collision energy, 35%; activation Q, 0.250; activation time, 30 ms) in 
the LTQ instrument after the accumulation of 10,000 ions. Maximal fill- 
ing times were 1,000 ms for the full scans and 150 ms for the MS/MS scans. 
Precursor ion charge state screening was enabled, and all unassigned 
charge states as well as singly charged species were rejected. The dynamic 
exclusion list was restricted to a maximum of 500 entries with a maximum 
retention period of 90 s and a relative mass window of 10 ppm. The lock 
mass option was enabled for survey scans to improve mass accuracy. The 
data were acquired using Xcalibur software. 

Peptide quantification. Quantification was performed with the Max- 
Quant program, version 1.0.7.4 (47), and was based on a two-dimensional 
centroid of the isotope clusters within each SILAC pair. The generation of 
the peak list, SILAC- and extracted ion current-based quantitation, calcu- 
lation of the posterior error probability, calculation of a false discovery 
rate based on search engine results, peptide-to-protein group assembly, 
and data filtration and presentation were carried out using the MaxQuant 
program. The derived peak list was searched with the Mascot search en- 
gine (version 2.1.04; Matrix Science, London, United Kingdom) against a 
concatenated database combining 80,412 proteins from the International 
Protein Index (IPI) human protein database, version 3.6 (forward data- 
base), and the reversed sequences of all proteins (reverse database). 

Data deposition. LC-MS/MS data from this project for the IBV N 
protein were deposited in the Proteomics Identifications (PRIDE) data- 
base using the PRIDE converter tool (48, 49) and also the IntAct database 
(50, 51) and were deposited through the IntAct database to the IMEx data 
resource and assigned the identifier IM-15828. 

Bioinformatic analysis. Ingenuity pathway analysis (IPA; Ingenuity 
Systems) and STRING analysis were used to analyze the cellular protein 
data sets and to group proteins into similar functional classes and examine 
potential interactions. IPA networks were generated using data sets con- 
taining gene identifiers and corresponding expression values, which were 
uploaded into the application. Each gene identifier was mapped to its 
corresponding gene object in the Ingenuity Pathways Knowledge Base 
(IPKB). Focus genes were overlaid onto a global molecular network de- 
veloped from information contained in the IPKB. Networks of these focus 
genes were then algorithmically generated on the basis of their connectiv- 
ity. Graphical representations of the molecular relationships between 
genes/gene products were generated. Genes or gene products are repre- 
sented as nodes, and the biological relationship between two nodes is 
represented as an edge (line). All edges are supported by at least one 
reference from the literature or from canonical information stored in the 
IPKB. Human, mouse, and rat orthologs of a gene are stored as separate 
objects in the IPKB but are represented as a single node in the network. 
STRING analysis was performed using version 9 (52), and STRING has 
been found to be a robust noncommercial algorithm for predicting pro- 
tein-protein interactions (53). 

Immunoblot analysis. Polyvinylidene difluoride (PVDF) membranes 
(Millipore) were activated in 100% methanol and equilibrated by immer- 
sion in SDS-PAGE transfer buffer (25 mM Tris-HCl, 192 mM glycine, 
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FIG 1 (A) Strategy for analyzing cellular interacting partners of IBV N protein via immunoprecipitation of an EGFP-N protein expressed in cells labeled by 
SILAC. Immunoprecipitations (IPs) were carried out using GFP-trap beads to minimize nonspecific binding. Labeling by SILAC was employed, as it allowed less 
stringent wash conditions to be employed, as contaminating proteins should be present in roughly equal amounts in both samples. Proteins showing an increased 
abundance with EGFP-N as the bait versus EGFP as the control bait are more likely to represent true interactions. Various stages in this process are highlighted. 
(B) Silver-stained SDS-polyacrylamide gel showing EGFP and EGFP-N pulldowns at various salt concentrations (indicated). (C). Immunoblot analysis of EGFP 
and EGFP-N pulldowns at various salt concentrations (indicated) to test the effects of buffer conditions on protein binding. Numbers to the left of the gels are 


molecular masses (in kilodaltons). 


20% [vol/vol] methanol). Transfer between the SDS-polyacrylamide gel 
and the membrane was performed using a Bio-Rad semidry blotting sys- 
tem according to the manufacturers’ instructions. Transfers were per- 
formed at 15 V for 60 min. After transfer, PVDF membranes were blocked 
in 10% (wt/vol) skimmed milk powder (Marvel) made up in Tris-buff- 
ered saline containing Tween 20 (TBS-T; 50 mM Tris-HCl, pH 8.3, 150 
mM NaCl, 0.5% [vol/vol] Tween 20) and incubated for 60 min. Blocking 
reagent was discarded, and primary antibody in TBS-T containing 5% 
(wt/vol) skimmed milk powder was added to the membrane and agitated 
for 60 min. Unbound antibody was removed by washing in TBS-T. Horse- 
radish peroxidase-conjugated secondary antibody in TBS-T containing 
5% (wt/vol) skimmed milk powder was then added to the membrane and 
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the membrane was agitated for 60 min. Unbound antibody was removed 
by washing in PBS-Tween 20. Protein-antibody complexes were visual- 
ized using an enhanced chemiluminescence (ECL) system (Amersham 
Biosciences). 

Ribosome cosedimentation. Ribosome cosedimentation was per- 
formed on 293T cells. At 2 h prior to harvest, cells were treated with 100 
g/ml cycloheximide. A total of 5 X 10° cells were used for each experi- 
ment. These cells were resuspended in 1 ml of cytoplasmic lysis buffer (10 
mM Tris-HCl [pH 7.4]. 2.5 mM MgCL, 100 mM NaCl, 0.5% [vol/vol] 
NP-40 alternative) and incubated on ice for 5 min before passing the 
sample through a 25-gauge needle three times. The samples were then 
centrifuged at 1,500 X g for 2 min and then 1,500  g for 15 min, with the 
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TABLE 1 SILAC ratios of proteins identified from the IBV N protein TABLE 1 (Continued) 
immunoprecipitation experiment* Binding No. of Sequetive 
Binding No. of Sequence Gene name ratio peptides coverage 
Gene name ratio peptides coverage DDX21 5.48 B 11 
RPL19 a8 3 18.4 RPL9 5.42 5 28.6 
RPL35 4.98 5 30.1 NONO 5.42 12 25.3 
CAPRINI1 4.32 9 12.8 SFRS1 5.34 6 29.8 
G3BP-2 4.02 11 29.5 RPL27 5.19 2 22.9 
RPL31 3.99 3 ey ie) SFPQ a9 10 17.5 
RPL13 3.45 2 9 HNRNPD 2.15 9 29.9 
RPS18 3.20 16 69.7 RPLIOA a.15 7 36.5 
RPS3A 2,23 14 48.9 RPL23A 5.00 3 17 
G3BP-1 2.21 14 38.8 DDX1 4.99 13 24.5 
ALYREF 1.94 2 6.8 XRCC5 4.97 12 19.7 
RPS9 1.89 3 95. HNRNPA/B 4.95 8 31.6 
RPS19 1.55 14 62.1 SSB 4.89 vi 20.1 
RPS16 0.67 i3 63 RPL5 4.79 6 26.6 
RPS5 0.26 1 42.6 HNRNPAI 4.69 16 43 
RPL24 9.81 4 21.4 GSK3B 4.69 5) 15.9 
RPS8& 9.76 8 35.9 C22orf28 4.60 11 27.5 
RPL26 9.59 8 48.4 RPLO 4.53 vs 28.4 
RPL22 9.54 4 48.4 RPL14 4Al 3 15 
YBX1 9.47 9 43 SFRS2 4.37 ) 10.9 
RPS11 9.45 10 54.4 HNRNPQ 4.35 13 26.5 
RPS10 9.39 8 33.9 CSDA 4.34 6 25.8 
RPS17 9.07 7 65.2 RPL30 4.34 6 58.3 
RPS25 9.07 8 46.4 DDX3X 4,32 10 18.3 
SFRS4 8.94 Z B22 RPL11 4.31 6 33.1 
RPS3 8.92 18 69.1 DDX17 4.31 16 22.8 
RPL4 8.74 7 19.7 DHX9 4.13 32 29.8 
RPS14 8.54 yi 40.4 FAM98A 4.10 2 12.1 
RPS26 8.40 4 37.4 HNRNPA3 4.07 6 25.7 
SFRS3 8.30 3 24.4 PABP1 4.05 34 48.1 
RPS23 8.15 5 29.4 DHX30 4.01 13 12.8 
RPS12 8.02 8 63.6 HNRNPU 3.96 21 25.2 
RPL15 8.00 4 21.6 FUS 3.92 4 10.8 
RPS13 7.97 8 45 SNRPD3 3.85 3 22.9 
HIST1H1D 719 > 17.2 RPL12 3.84 6 48.5 
RPS29 72 5 29.6 SERBP1 3.82 10 33.1 
GSK3A 7.60 5. 13 BLM ae 2 1.6 
RPSI15A 7.60 8 46.9 RPSA 377 10 38 
RPL13A 7.45 3 16.7 HNRNPR a7 5 9 
RPS15 7.43 4 22.1 HNRNPA2/B1 3.75 1S 45 
RPL23 7.24 3 27.1 NPMI1 3.07 10 371 
RPL22L1 7.19 ) 51.6 RPLP1 3.56 2 28.9 
RPL7 7.03 v 25.7 LARPI 3.45 28 31.5 
RPLI8A 6.96 4 22.2 RBM14 3.40 ] 3.6 
RPL10 6.94 7 35.2 SFRS7 3.09 2 8.8 
RPS7 6.88 8 47.4 YTHDF2 3.38 i} 10.2 
RPLI8A 6.78 4 27.8 GNB2L1 Ber) 13 50.1 
RPS4 6.61 14 48.3 HNRNPG 3.30 2 6.6 
NCL 6.58 22 29.7 IGF2BP1 3.29 20 40.2 
RPL21 6.38 2 11.2 PTBP1 3.27 Lt 26.8 
RPL17 6.33 5 32.1 PABP4 320 30 41.5 
RPS24 6.33 2 20.3 DHX36 and 6 7 
DDX5 6.29 14 23.8 MOV10 3.21 3 3.6 
RPL8& 6.26 3 13.4 DHX57 3.18 2 1.7 
HINRNPAO 6.26 4 17.7 HIFX 3.11 3 21.1 
cl4orf166 5.91 3 11.9 ATXN2L 3.09 3 2.8 
RPL3 5.83 4 8.9 FMRI 3.05 3 4.9 
RBM3 5.76 2 19.1 DHX15 3.01 2 2.9 
USP10 5.55 5 9.1 HNRNPHI1 3.01 6 19.5 
RPL27A S51 3 aL HINRNPF 2.98 3 10.6 
4 37.9 


XRCC6 5.51 15 28.1 H2AFX 2.97 


(Continued on following page) 


September 2013 Volume 87 Number 17 jviasm.org 9489 


Aresqr] [ed1Ipey SION OSN AQ L102 ‘v [dy uo /Bio-wseIAl//:dyjy WO. pepeojumMoq 


Emmott et al. 


TABLE 1 (Continued) 


Binding No. of Sequence 

Gene name ratio peptides coverage 
HNRNPK 2.95 12 35.6 
ILF2 2.91 3 10.3 
DDX6 2.80 2 4.3 
LRRC59 2.78 3 11.4 
HNRNPM 2.69 10 17.3. 
HDLBP 2.68 9 8.9 
HNRNPH3 2.54 2 8.4 
PARP1 2.50 7 8.6 
ELAV1 2.36 5 16.4 
HNRNPC 2.36 5 19.6 
ILF3 2.35 13 16.3 
MYL6 2:59 2 16.6 
RPS27A 2.30 > 42.3 
ZCCHC3 2329 2 6.4 
HNRPDL 2.24 5 11.7 
LRPPRC 2.15 12 11 
MYH9 2.13 34 22.5 
IGF2BP3 2.11 12 29) 
HNRNPL 2.07 2: 4.1 
SND1 2.05 2 2.5 
RENTI1 2.03 20 22 


* Ratios are organized from the potential strongest interaction to the weakest 
interaction. Shown are the gene name of the protein, binding ratio (EGFP-N/EGFP), 
number of different peptides used to identify the protein, and, from this value, the 
sequence coverage (which is calculated by dividing the number of amino acids observed 
by LC-MS/MS by the protein amino acid length). Proteins identified to have one 
peptide and/or a binding ratio of less than 2.00 have been excluded. 


supernatant retained at each stage. Sucrose gradients (5% and 20%) were 
made up in 10 mM Tris-HCl, pH 7.4, 100 mM NaCl, and 20 mM EDTA 
and passed through a 22-um-pore-size filter. These were then used to 
generate 12-ml 5 to 20% sucrose gradients using a Hoefer 15 gradient 
maker. The cytoplasmic fraction (200 wl) was then layered over the gra- 
dient and centrifuged at 287,000 X gin a swinging-bucket rotor for 5 h at 
4°C. Fractions (1 ml) were collected from the bottom of the tube by needle 
puncture. 

Virus infection. Cells of the Vero cell line (an African green monkey 
kidney-derived epithelial cell line) were maintained in DMEM supple- 
mented with 10% fetal calf serum (FCS). IBV Beaudette US, a Vero cell- 
adapted isolate of Beaudette CK (28), was propagated in Vero cells. All cell 
culture experiments in this study were conducted in the absence of anti- 
biotic or antifungal agents and were performed on actively replicating 


i) 5 10 61S 


subconfluent cells (i.e., cells not undergoing contact inhibition). At 60% 
confluence, cells were infected with IBV and then incubated for 1 h at 
37°C, after which the cells were incubated in maintenance medium; cells 
were then processed for analysis at 12 h postinfection. 

Cell viability assay. Cell viability/cytotoxicity during siRNA treat- 
ment was assessed using a colorimetric 3-(4,5-dimethylthiazol-2-yl)-2,5- 
diphenyltetrazolium bromide (MTT) assay. Vero cells were seeded at a 
density of 3.8 X 10° (for 40 to 50% confluence) cells/well in clear 96-well 
microplates at 24 h prior to small interfering RNA (siRNA) treatment, 
following which 24 mg MTT powder was dissolved in 10 ml DMEM (plus 
10% [vol/vol] FBS) at 37°C to make a 10 mM solution and the solution 
was filtered through a 0.2-m-pore-size filter. The medium was aspirated, 
and the wells were washed with PBS. The absorbance at between 425 and 
570 nm was then measured on a Dynex plate reader. A dose-response 
curve was then produced to assess cell viability. 

siRNA-mediated knockdown. For siRNA-mediated knockdown of 
target mRNAs, the following approach was used: two sequences for the 
same target gene were used alongside a negative control (Qiagen). The 
lyophilized oligomers were resuspended in nuclease-free water to make a 
10 uM stock solution, and aliquots were further diluted to provide 10 nM 
working stocks. Lipofectamine 2000 was used (according to the manufac- 
turer’s instructions) to transfect Vero cells with siRNA. For a typical 24- 
well plate transfection, cells were seeded at a density of 3 X 10* in antibi- 
otic-free DMEM at 24 h prior to transfection to achieve 30 to 40% 
confluence. Prior to transfection, the growth medium was replaced with 
400 yl fresh antibiotic-free DMEM and 10 pmol of each siRNA pair (20 
pmol for the negative control) and the DMEM and siRNA were mixed 
with 50 pl well Opti- MEM medium (Invitrogen) at room temperature for 
5 min. In parallel, 1 jl/well Lipofectamine 2000 was mixed with 50 pl/well 
Opti-MEM medium at room temperature for 5 min. The dilute oligomers 
and the dilute Lipofectamine 2000 were then mixed gently and incubated 
at room temperature for 20 min. The oligomer-Lipofectamine 2000 com- 
plexes were then added to each well (100 «l/well) dropwise to give a final 
siRNA concentration of 20 pmol total siRNA/100 pl Opti- MEM medium. 
Control wells of reagent only and medium only were included to assess 
any background reactivity. 

Confocal imaging. Cells were fixed using a 4% (wt/vol) solution of 
paraformaldehyde made up in PBS and incubated for 10 min at 20°C. 
Following fixation, cells were made permeable with 0.2% (vol/vol) Triton 
X-100 made up in PBS for 10 min at 20°C. Primary antibody was added in 
PBS containing 1% FBS, and primary and secondary antibodies were 
tested to ensure no cross-reactivity. Slides were imaged on an LSM 510 
META confocal microscope (Carl Zeiss). All images were captured using 
either a X40 or a X63 objective and a digital zoom factor of 1 to 4 within 
the software. The rainbow feature of the software was used to ensure that 
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FIG 3 STRING (version 9.0) analysis of the cellular proteins identified by the SILAC pulldown approach to be potentially interacting with IBV N protein. 
Proteins group into two large units: those associated with translation (right) and those involved in RNA modification/processing (left). 


images were taken wherever possible within the linear range (unless stated 
otherwise). 


RESULTS 


To investigate the cellular interactions of the IBV N protein, a 
SILAC-based quantitative proteomic approach coupled to an im- 
munoprecipitation strategy based on expression of an EGFP-N 
fusion protein in cells was used (Fig. 1A). Expressing the target 
protein as an EGFP-fusion protein has been shown to improve 
sensitivity and allow discrimination of specific from nonspecific 
interactions with the target protein (54, 55). 

Optimization of immunoprecipitation conditions. From 
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previous experience investigating proteome changes in IBV-in- 
fected cells using various model cell lines (Vero [primate-derived] 
and DF1 [avian-derived] cells) (45, 46), a human cell line was used 
in this study to provide cellular proteins with potential interaction 
with N protein. This was due to the superior annotation of the 
human genome allowing protein identification. 293T cells were 
also chosen because of their efficiency in taking up exogenous 
plasmid DNA by calcium phosphate (>45% by fluorescence-ac- 
tivated cell sorting; data not shown) compared to that of Vero or 
DFI1 cells, thus allowing higher exogenous protein expression, 
which could be monitored using the EGFP tag. IBV N protein was 
expressed in 293T cells as part of an EGFP fusion protein (EGFP- 
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N), as previously described (41-44), and the GFP-trap system was 
then applied to allow isolation of potential interaction partners 
from the cell lysates. 

Various concentrations of NaCl (150, 300, and 500 mM) were 
tested in the wash buffer to determine the optimum binding con- 
ditions. Immunoprecipitated samples were analyzed by SDS- 
PAGE, and proteins were visualized by silver staining (Fig. 1B). 
EGFP alone gave heavy staining between the 28- and 38-kDa 
marker bands, with a few faint bands being present at different 
molecular masses under the 150 mM NaCl wash condition. These 
products became less apparent as the concentration of NaCl was 
increased. In contrast, immunoprecipitation from the EGFP-N 
samples indicated a large number of heavily stained products from 
alow molecular mass (14 kDa) upwards. A densely staining band 
at approximately 75 kDa was identified as the EGFP-N fusion 
protein. At either a 150 mM or a 300 mM NaCl concentration, 
silver-stained bands in the EGFP-N samples did not show signif- 
icant decreases in their intensity. However, at 500 mM NaCl, in- 
tensity differences became more apparent. Immunoblot analysis 
using either an anti-EGFP antibody or an anti-IBV antibody con- 
firmed the expression and immunoprecipitation of EGFP and 
EGFP-N, respectively (Fig. 1C). Note that the species visualized 
migrating between the 98- and 64-kDa markers in the anti-GFP 
and anti-IBV EGFP samples was present in other pulldown sam- 
ples, and we postulate that it represents nonspecific binding of a 
host protein by the GFP-trap, with the apparent high level of en- 
richment leading to antibody cross-reactivity. Taken together, a 
300 mM NaCl concentration was used in the wash buffer, as this 
appeared to yield a decrease in nonspecific binding in the EGFP 
pulldown samples while maintaining staining patterns in the 
EGFP-N samples. Further, SILAC coupled to LC-MS/MS was 
used to discriminate between background and specific interac- 
tions with the N moiety. 

LC-MS/MS analysis. SILAC-labeled samples (ROKO EGFP, 
R6K4 EGFP-N) were immunoprecipitated, and equal volumes 
were combined and analyzed by LC-MS/MS to identify and quan- 
tify target proteins (45, 56, 57). Out of the 245 cellular proteins 
identified, 232 were also quantified. This data set and the associ- 
ated mass spectra were uploaded to the PRIDE repository (49). 
For many quantitative proteomics analyses, an arbitrary 2-fold 
cutoff has been used to select potential true versus nonspecific 
interactions. When applied to this data set, 142 cellular proteins 
met this criterion (Table 1). These proteins were used in down- 
stream analysis and also deposited with the IntAct database. 

Bioinformatic analysis of the interactome data set. Ingenuity 
pathway analysis was used to investigate whether the identified 
proteins grouped into distinct functional clusters. This revealed 
that the identified proteins were mainly involved in protein syn- 
thesis and RNA posttranscriptional modification (Fig. 2). Inde- 
pendent analysis of the potential network interactions using 
STRING revealed clustering of proteins into those involved in 
translation (Fig. 3, right) and RNA processing/modification (Fig. 
3, left). Many proteins were associated with the large (28 proteins) 
and small (21 proteins) ribosomal subunits, as well as 16 hnRNP 
proteins with other groups, including the DDX and DHX RNA 
helicases and splicing factors and the kinases GSK3a and GSK3B. 

The efficiency of binding to EGFP-N over that to EGFP does 
not reflect the abundance of the proteins in the cellular pro- 
teome. One hypothesis that could explain the data set was that the 
degree of enrichment of binding to EGFP-N merely reflected 
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FIG 4 Immunoblotting confirmation of the interactions of identified cellular 
proteins with EGFP-N protein. (A) Agarose electrophoresis analysis of a small 
aliquot of the input lysate (to be used in the pulldown assay) in the absence (—) 
and presence (+) of RNase. Lane M, a 100-bp DNA ladder. (B) The pulldown 
experiment was repeated in the presence and absence of RNase and immuno- 
blot analysis of selected cellular proteins. 


the abundance of the specific protein in the cellular proteome; i-e., 
the most abundant cellular proteins would be overrepresented. To 
investigate this, we compared the EGFP-N data set with data in the 
PaxDb: Protein Abundance across Organisms database, which 
ranks proteins according to their abundance in various human cell 
types on the basis of spectral counting data. The eight most abun- 
dant proteins in human cells at the time of analysis were recorded 
as APOA2, RBP4, APOC2, ALB, TTR, APOAI, APOC1, and 
APOC3; however, none of these were part of the 142 proteins 
associated with EGFP-N. RPL19 was the most enriched protein 
out of the 142, and this is ranked the 446th most abundant out of 
12,797 cellular proteins. RENT1, which was the least enriched 
protein out of the 142 selected, is ranked 1,564 out of 12,797. 
Between these two proteins were proteins that occurred with dif- 
ferent frequencies. For example, GSK3a and GSK3B are ranked 
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FIG 5 Cosedimentation of EGFP-N protein with the small ribosomal subunit. Cosedimentation was performed on a 5 to 20% sucrose gradient in the presence 
of EDTA to induce separation of the large and small ribosomal subunits. The migration of these subunits was determined by extraction of RNA from the fractions 
and running of these on a 1% agarose gel containing ethidium bromide (EtBr), allowing visualization of the 18S and 28S rRNAs. Migration of the EGFP control 
or EGFP-N protein was determined by immunoblotting using an anti-GFP antibody. 


the 3,755th and 2,962nd most abundant out of 12,797 proteins, 
respectively, and G3BP-2 is ranked 1,568 out of 12,797. Therefore, 
there is no apparent correlation between the abundance of a pro- 
tein in a cell and its association with N protein. 

Validation of SILAC immunoprecipitation results by immu- 
noblotting. To increase confidence in the SILAC EGFP-N pull- 
down results, the experiment was repeated in the absence of label 
and immunoblotting was used to confirm the identity of N-inter- 
acting cellular proteins using specific antibodies (Fig. 4). As N 
protein is a known RNA-binding protein (24), analysis of the 
EGEP-N protein input lysate (to be used in the pulldown assays) 
indicated the presence of RNA (Fig. 4A), and this interaction may 
also determine an association with some cellular proteins. Exam- 
ples of proteins used for validation from the data set were taken 
from Table 1. All of the proteins were confirmed to be interacting 
with the EGFP-N protein, with the interaction with NONO and 
poly(A)-binding protein (PABP) being potentially mediated by 
RNA (Fig. 4B). Other proteins, such as nucleolin (NCL), have 
previously been shown to interact with native IBV N protein in an 
RNA-independent manner (40). No identified proteins were 
shown to interact with EGFP (Fig. 3B). There is also some evi- 
dence from the EGFP blot that some EGFP may be cleaved from 
EGEP-N, and this may be through the action of caspase (58). 

Interaction of IBV N with the small ribosomal subunit. As 
over a third of the identified cellular proteins were ribosomal pro- 
teins, a potential interaction between N protein and ribosomal 
proteins was suggested. To investigate whether EGFP-N protein 
sedimented with either or both of the ribosomal subunits, 28S and 
18S rRNAs were used to show the migration of the large (60S) and 
small (40S) ribosomal subunits, respectively (Fig. 5). In EGFP- 
expressing control cells, EGFP sedimented at the top of the gradi- 
ent in fractions 9 to 12 and was most concentrated in fraction 11, 
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while the ribosomal subunits occupied fractions 2 to 7. In con- 
trast, EGFP-N protein was the most concentrated in fraction 7, the 
least concentrated in fraction 8, and at an intermediate level in 
fraction 6, a pattern duplicated by the 185 rRNA. This suggested a 
preference for one of the ribosomal subunits. The small ribosomal 
subunit is associated with translation initiation. 

Interaction of the EGFP-N protein with G3BP. Three pro- 
teins with some of the highest potential binding ratios were 
caprin-1, G3BP-1, and G3BP-2 (Table 1) and are involved in the 
formation of cytoplasmic stress granules (59, 60). Indirect immu- 
nofluorescence confocal microscopy was used to investigate 
whether the EGFP-N protein colocalized with G3BP. This pro- 
vided an independent validation of the interaction between these 
proteins separate from immunoblotting. To induce stress granule 
formation and, hence, G3BP, cells were treated with sodium ar- 
senite (Fig. 6) and G3BP was used as a stress granule marker. In 
cells lacking stress granules, G3BP staining was diffuse and cyto- 
plasmic, whereas in arsenite-stressed cells, G3BP was concen- 
trated in punctate cytoplasmic foci. In EGFP-expressing cells, 
arsenite treatment successfully induced stress granule formation, 
with granules showing no enrichment of EGFP. In contrast, in 
EGFP-N protein-expressing cells, arsenite treatment caused 
EGFP-N protein to relocalize to G3BP-rich foci (Fig. 6, yellow 
signal), suggestive of an interaction between this protein and 
EGFP-N protein. 

Potential relocalization and association of cellular proteins 
with N protein in IBV-infected cells. To investigate whether any 
of the cellular proteins that interacted with N protein were either 
redistributed in IBV-infected cells or colocalized with N protein, 
Vero cells were infected with IBV and subcellular localization was 
determined at 12 h postinfection (Fig. 7). Selection of cellular 
proteins was on the basis of a potential association and also the 
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FIG 6 Indirect immunofluorescent confocal microscopy showing colocalization of EGFP-N protein with the stress granule marker in G3BP in arsenite-stressed 
293T cells. Cells were mock treated (—) or treated (+) with 0.5 mM sodium arsenite for 1 h prior to fixation. Nuclei are colored blue (DAPI [4’,6-diamidino- 
2-phenylindole]), EGFP is in green, G3BP is in red, and a merge image is presented. Colocalization is suggested by a yellow signal. Bars, 10 mm. 


availability and avidity of antibodies. The ribosomal protein L19 
(shown in red) localized to both the cytoplasm and the nucleolus 
both in bystander cells and also in infected cells, as would be pre- 
dicted for a ribosomal protein. Additionally, the data indicated 
that L19 colocalized with N protein in infected cells (green), as 
shown by the yellow signal in the merged image. PABP (red) co- 
localized with N protein (green) in infected cells (yellow in the 
merge image). Both B23.1 (a nucleolar protein, indicated in red) 
and hnRNPA2/B1 (a nuclear protein, indicated in red) localized to 
the nucleolus and nucleus, respectively, in bystander cells. How- 
ever, in infected cells (green), B23.1 (red) and hnRNPA2/B1 (red) 
localized not only to the nucleolus and nucleus but also to the 
cytoplasm, as shown by the red signal in the cytoplasm. Colocal- 
ization was not possible to investigate, given the difference in sig- 
nal intensity in the cytoplasm between B23.1/hnRNPA2/B1 and N 
protein. In contrast, no relocalization or colocalization with N 
protein was observed with hnRNPA1, despite capturing the signal 
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for hnRNPA1 beyond the linear range (i.e., red-saturated pixels). 
Likewise, no relocalization or colocalization with N protein was 
observed with caprin-1, NONO, or ALY in IBV-infected cells. 
Potential importance of the cellular proteins that associate 
with N protein in the viral life cycle. Several approaches have 
been used to investigate the importance of cellular proteins in a 
virus life cycle that interact with viral proteins; generally, these 
involved ablation or partial ablation of the cellular protein using 
siRNA (e.g., for IBV [61]). To investigate whether the identified 
cellular proteins that interacted with EGFP-N were involved in the 
replication of IBV, Vero cells were transfected with two specific 
siRNAs per mRNA directed toward selected cellular targets, prior 
to infection with IBV. The cellular targets were chosen on the basis 
of demonstrated siRNA knockdown and potential association 
with the N protein, including proteins of high and low probability. 
Four target proteins that matched these criteria were chosen: 
RPL19, GSK3a/GSKB, nucleolin, and procyclic acidic repetitive 
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FIG 7 Indirect immunofluorescent confocal microscopy analysis of the localization of selected cellular proteins that potentially interact with N protein in 
IBV-infected cells. Cellular proteins are shown in red, IBV-infected cells are shown in green, and a merged image is presented. Colocalization is suggested by a 


yellow signal. Bars, 10 mm. 


protein 1 (PARP1). These proteins had EGFP-to-EGFP-N bind- 
ing ratios of 15.78, 7.60/4.69, 6.58, and 2.5, respectively. In addi- 
tion, cells were transfected with a nontargeting siRNA pair. West- 
ern blotting indicated that the abundance of all four proteins was 
decreased in the presence of the specific siRNA but not the non- 
targeting siRNA (Fig. 8A). However, the efficiency of ablation was 
different, with RPL19, PARP1, and GSK3 being the most effi- 
ciently depleted. Western blotting indicated that in cells depleted 
of nucleolin, RPL19, GSK3a, and GSK38, the abundance of N 
protein decreased compared to that for the control (e.g., nontar- 
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geting siRNA-treated cells) or untreated cells (Fig. 8A). The ob- 
served decreased abundance of N protein was not due to any gross 
effects on cell viability (Fig. 8B), although there was an approxi- 
mately 15 to 20% greater reduction in cell viability in GSK3-de- 
pleted cells than in the other treated cells. Given the decrease in the 
abundance of N protein in infected cells in which the abundance 
of nucleolin, RPL19, and GSK3 was decreased, we hypothesized 
that synthesis of viral RNA may be affected. To test this hypothe- 
sis, quantitative reverse transcription-PCR (qRT-PCR) using a 
specific primer pair to identify the genomic RNA was used to 
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investigate the abundance of this molecule under the different 
conditions (Fig. 8C). The data indicated that the amount of 
genomic RNA was decreased in infected cells depleted of nucleo- 
lin, RPL19, or GSK3 but not in infected cells depleted of PARP1 or 
in untreated cells (Fig. 8C). 


DISCUSSION 

This study represents the first work on the use of SILAC-based 
LC-MS/MS to study the interaction of the coronavirus N protein 
with host cell proteins and to tentatively identify the cellular in- 
teractome of this protein. The SILAC-based LC-MS/MS approach 
allowed the discrimination between cellular proteins that bound 
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both to EGFP and/or the binding matrix and to N protein. Prior to 
this work, the sole known interactions confirmed for the IBV N 
protein were an interaction with the nucleolar protein nucleolin 
and colocalization (but not necessarily interaction) with fibrillarin 
(40). Studies carried out with different coronavirus N proteins 
demonstrated interactions with the kinase GSK3B (62), as well as 
the A-type hnRNPs (63), which offered the potential for conser- 
vation across the Coronaviridae. The approach described herein 
identified a potential of 142 cellular partners, including a large 
number of hnRNPs, ribosomal proteins, and the kinases GSK3a 
and GSK38, and confirmation of the interaction with nucleolin as 
well as a novel interaction with B23.1 and with several stress gran- 
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FIG 8 Investigation of the proviral/antiviral activity of selected cellular proteins that potentially interact with N protein in IBV-infected Vero cells. (A) 
Representative immunoblots showing the abundance of the target protein and IBV N protein in mock-infected and infected Vero cells either treated with siRNAs 
targeting the specific mRNA of interest or treated with a nontarget siRNA control. (B) MTT cell viability assay of Vero cells treated with either the nontarget 
siRNA or the specific siRNA to the mRNAs encoding the selected protein of interest. The experiment was performed three times in triplicate for each condition. 
(C) qRT-PCR analysis of the abundance of the IBV genomic RNA in either mock-infected cells (mock) or cells treated with the nontarget siRNA or the siRNA 
specific to the mRNAs encoding the selected protein of interest. Replication efficiency is related to the level observed in untreated cells. si-non target, nontargeting 
siRNA; si-Nucleolin, siRNA targeting nucleolin; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; si-PARP1, siRNA targeting PARP1. 


ule marker proteins known to form a heterodimeric complex. N 
protein is unlikely to interact individually with all of the identified 
proteins, and these data are most likely to represent those for 
many complexes of cellular proteins that bind via a number of hub 
proteins, such as nucleolin (35). 

Comparison to other viral proteins where we have used this 
approach, including the human respiratory syncytial virus 
(HRSV) NS1 and N32 proteins (64), indicated that the IBV N data 
set was unique in terms of highlighting specific proteins and 
classes of proteins that potentially bound to the N protein. For 
example, ribosomal proteins were not the top proteins with bind- 
ing to either HRSV NS1 or NS2. Similar to the analysis of the NS1 
cellular interactome (64), the data indicated that the SILAC-based 
LC-MS/MS approach coupled to the GFP-trap immunoprecipita- 
tions could be used to identify protein complexes that bound to 
the viral protein of interest. Analysis of the related arterivirus 
PRRSV N protein (65, 79) indicated interaction with translation 
initiation factors, proteins involved in mRNA stability, hnRNPs, 
and DEAD RNA helicases. Although PRRSV N protein interacted 
with some ribosomal proteins, this was not necessarily the reper- 
toire of ribosomal proteins associated with IBV N protein. 

Results of repeat pulldown assays with EGFP-N protein con- 
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firmed the data in the mass spectrometry-based data set in the 
absence of label using immunoblotting and also confocal micros- 
copy. Gene depletion analysis on selected targets using siRNA was 
used to demonstrate the potential functional relevance of cellular 
proteins in the biology of the virus, although off-target effects 
could not be ruled out. For example, depletion of GSK-3 in IBV- 
infected cells resulted in a decrease in the abundance of N protein 
and a reduction in the amount of genomic RNA compared to the 
results for the control treated cells. Although cells were generally 
viable, GSK3a and GSK3B are major cellular kinases and deple- 
tion of these molecules also likely affected other cellular pathways. 
However, GSK3 has previously been shown to be responsible for 
phosphorylation of the SARS-CoV and MHV N proteins (62). In 
addition, treatment of IBV-infected cells with lithium chloride 
(LiCl), which is a GSK3 inhibitor, resulted in decreased amounts 
of viral RNA and reduced viral titers (66, 67). 

Interactions between the coronavirus N protein and hnRNPs 
have previously been reported for other members of the Corona- 
viridae, including SARS-CoV and MHV N proteins (68, 69), and 
have been shown through overexpression/mRNA depletion anal- 
ysis to potentially be involved in coronavirus replication and tran- 
scription (70-72). The discovery that these hnRNP molecules 
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bind N protein in an RNA-independent manner lends credence to 
the hypothesis that these proteins may be involved in replication 
and transcription of viral RNA, with at least some of the molecules 
able to functionally substitute for others in the event of their loss 
(72). 

Ribosomal proteins also represented a large proportion of the 
cellular proteins identified in this analysis. Very limited data on 
the role of N protein in translation are available, but in MHV, N 
protein was found to favor the translation of viral over cellular 
transcripts, in a mechanism requiring N protein to bind to telo- 
mere repeat sequences (73, 74). There was a 4-fold increase in the 
efficiency of translation for reporter sequences containing MHV 
leader RNA compared to that for control sequences containing 
human alpha-globin or reverse complement leader sequence. 
EGEP-N was also shown to interact with various RNA helicases. 
One of these, DDX1, has also been shown to interact with IBV and 
SARS-CoV nonstructural protein 14 (61). Depletion of DDX1 in 
IBV-infected cells resulted in a reduction of IBV replication (61). 

Several of the cellular proteins identified as interacting with the 
IBV N protein have also shown to be essential cofactors in the 
replication of other RNA viruses, indicating the importance of 
recruitment of cellular proteins in the life cycle of viruses. For 
example, Cl4orf166, which is involved in the modulation of 
mRNA transcription by polymerase I, has been shown to interact 
with the PA subunit of the influenza virus polymerase, and silenc- 
ing causes a reduction in polymerase activity and a reduction in 
virus titer (75). DDX17, DDX5, NPM1, and HNRNPM were 
shown to be required for efficient activity of the influenza virus 
polymerase (76). The rubella virus capsid proteins binds to PABP 
and inhibits cellular protein synthesis (77). CAPRIN1 associates 
with the Japanese encephalitis virus core protein, and this inter- 
action facilitates virus growth as well inhibits stress granule for- 
mation (78). 

In summary, this study identified and characterized the poten- 
tial cellular interactome of the coronavirus IBV N protein, using a 
SILAC-based quantitative proteomic methodology to allow dis- 
crimination of specific versus nonspecific interactions during the 
pulldown process. This interactome both confirms previous ob- 
servations made with other coronavirus and IBV N proteins with 
both overexpressed proteins and infectious virus and also pro- 
vides novel data that can be exploited to understand the interac- 
tion between the virus and the host cell. 
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